Draft version January 20, 2013 

Preprint typeset using I^T^X style emulateapj v. 11/10/09 



TORUS AND AGN PROPERTIES OF NEARBY SEYFERT GALAXIES: RESULTS FROM FITTING IR 
SPECTRAL ENERGY DISTRIBUTIONS AND SPECTROSCOPY 

Almudena Alonso-Herrero^'^, Cristina Ramos Almeida^, Rachel Mason^, Andres Asensio Ramos^'^, Patrick F. 
RocHE^, Nancy A. Levenson^, Moshe Elitzur^, Christopher Packham^°, Jose Miguel Rodriguez Espinosa^'^, 
Stuart Young^\ Tanio Diaz-Santos^^, and Ana M. Perez Garcia^'^ 

Draft version January 20, 2013 

ABSTRACT 

We used the CLUMPY torus models and a Bayesian approach to fit the infrared spectral energy 
distributions (SEDs) and ground-based high- angular resolution mid-infrared spectroscopy of 13 nearby 
Seyfert galaxies. This allowed us to put tight constraints on torus model parameters such as the 
viewing angle the radial thickness of the torus F, the angular size of the cloud distribution atorus, 
and the average number of clouds along radial equatorial rays Nq. We found that the viewing angle 
i is not the only parameter controlling the classification of a galaxy into a type 1 or a type 2. In 
principle type 2s could be viewed at any viewing angle i as long as there is one cloud along the line 
of sight. A more relevant quantity for clumpy media is the probability for an AGN photon to escape 
unabsorbed. In our sample, type Is have relatively high escape probabilities, Pesc 12 — 44%, while 
type 2s, as expected, tend to have very low escape probabilities. Our fits also confirmed that the tori 
of Seyfert galaxies are compact with torus model radii in the range 1-6 pc. The scaling of the models to 
the data also provided the AGN bolometric luminosities Lboi(AGN), which were found to be in good 
agreement with estimates from the literature. When we combined our sample of Seyfert galaxies with 
a sample of PC quasars from the literature to span a range of Lboi(AGN) 10^^ — lO^^ergs"^, we 
found plausible evidence of the receding torus. That is, there is a tendency for the torus geometrical 
covering factor to be lower (/2 ~ 0.1 — 0.3) at high AGN luminosities than at low AGN luminosities 
(/2 ^ 0.9 — 1 at 10^^~^^ ergs~^). This is because at low AGN luminosities the tori appear to have 
wider angular sizes (larger atorus) and more clouds along radial equatorial rays. We cannot, however 
rule out the possibility that this is due to contamination by extended dust structures not associated 
with the dusty torus at low AGN luminosities, since most of these in our sample are hosted in highly 
inclined galaxies. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — galaxies: structure — infrared: galaxies 



1. INTRODUCTION 

The unified model for active galactic nuclei (AGN) 
proposes the ubiquitous presence of an obscuring torus 
around their nuclei, with type 1 and type 2 AGN being 
intrinsically similar (Antonucci 1993). The central region 
of an AGN (including the Broad Emission Line Region, 
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BLR) is obscured when viewed along directions close to 
the equatorial plane of the torus, and it is then classi- 
fied as a type 2 AGN. In type 1 AGN the viewing angle 
is close to the polar direction of the torus, and thus we 
have a direct view of the central engine. This model re- 
ceived strong support from the fact that broad lines have 
been revealed in the spectra of the polarized emission of 
a number of type 2 AGN (e.g., Antonucci & Miller 1985; 
Tran, Miller, & Kay 1992). 

There is also indirect evidence of the presence of the 
dusty torus advocated by the unified model. First is the 
similarity of Seyfert Is and Seyfert 2s with respect to 
a given isotropic indicator of the AGN luminosity such 
as, [Olll]A5007, infrared (IR), hard X-ray, and radio lu- 
minosities of Seyferts (see e.g., Mulchaey et al. 1994; 
Alonso-Herrero, Ward, & Kotilainen 1997; Nagar et al. 
1999). The UV and soft X-ray continua of Seyfert 2s, on 
the other hand are underluminous relative to the type 
Is, because they are not transmitted through the torus. 
Second, the presence of double-cone morphology of the 
narrow line region (NLR) structures (e.g., Pogge 1989; 
Wilson & Tsvetanov 1996; Schmitt et al. 2003), with 
orientations similar to those of radio jets (Nagar et al. 
1999 and references therein), are interpreted as emission 
collimated by the torus. 

The role of the extinction produced by the host galaxy 
cannot be understated (Alonso-Herrero et al. 2003). 
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For instance, for some Seyfert 2s broad lines are iden- 
tified in the near-IR (e.g., Blanco, Ward, & Wright 1990; 
Ruiz, Rieke, & Schmidt 1994; Goodrich et al. 1994; 
Veilleux, Goodrich, & Hill 1997; Ramos Almeida, Perez- 
Garcia, & Acosta-Pulido 2009), and counter-cones are 
seen in direct and/or polarized near-IR light (e.g.. Pack- 
ham et al. 1997 for NGC 1068; Maiolino et al. 2000 for 
Circinus). Finally, there is the tendency of intermediate 
types (Seyferts 1.8+1.9) to be hosted in edge-on galaxies 
(Maiolino & Rieke 1995). 

In recent years much progress has been made to- 
ward understanding the properties of the molecular dusty 
torus. From the theoretical point of view, torus models 
with smooth density distributions (e.g.. Pier & Krolik 
1993; Granato & Danese 1994; Efstathiou & Rowan- 
Robinson 1995) have been superseded by more com- 
plex clumpy dust distributions (e.g., Nenkova, Ivezic, & 
Elitzur 2002; Honig et al. 2006; Nenkova et al. 2008a, 
2008b; Schartmann et al. 2008; Honig & Kishimoto 
2010). The clumpy models reproduce well the near-IR 
and mid-IR emission of different types of AGN (e.g.. Ma- 
son et al. 2006, 2009; Nenkova et al. 2008b; Schartmann 
et al. 2008; Polletta et al. 2008; Mor, Netzer, & Elitzur 
2009; Thompson et al. 2009; Nikutta, Elitzur, & Lacy 
2009; Ramos Almeida et al. 2009 and 2011, RA09 and 
RAll hereafter; Honig et al. 2010; Lira et al. 2011), 
and overcome some of the difficulties faced by smooth 
torus models in fitting IR data of AGN (see e.g., Alonso- 
Herrero et al. 2001, 2003). 

From the observational point of view, mid-IR inter- 
ferometric observations have revealed that the torus is 
relatively compact, with typical sizes of a few par sees, ge- 
ometrically thick, and with evidence of dumpiness (e.g., 
Jaffe et al. 2004; Tristram et al. 2007, 2009; Burtscher et 
al. 2009; Raban et al. 2009). Also, direct mid-IR imag- 
ing observations of Circinus gave a torus size of < 4pc 
(Packham et al. 2005). In terms of their molecular gas, 
the tori contain ~ 10^ M©, extend for tens of parsecs 
(typically radii of 30 pc) and have gas column densities, 
as derived from near-IR molecular hydrogen lines, rang- 
ing from 1 to 10 X 10^^ cm~^ (Davies et al. 2006; Hicks 
et al. 2009). 

This is the third paper in a series using high angular 
resolution IR observations, the CLUMPY torus models 
of Nenkova et al. (2008a, b), and a Bayesian approach 
for fitting the data (Asensio Ramos & Ramos Almeida 
2009) to derive the torus and AGN properties. In the 
first two papers of the series (RA09 and RAll) we fitted 
the IR photometric spectral energy distributions (SEDs) 
of a sample of nearby Seyfert 1 and Seyfert 2 galaxies. 
We constrained several torus model parameters, namely 
its angular width and the average number of clouds along 
radial equatorial rays, as well as the viewing angle to the 
torus. We found, in clear contrast with the simplest uni- 
fied model predictions, that type 2s appear to have tori 
with wider cloud distributions (the torus angular width 
c^torus parameter, see Figure 1, and Section 3.1) and more 
clumps than those of type Is. This may suggest that 
some of the properties of the tori of type 1 and type 2 
AGNs are intrinsically different. For a sample of nearby 
active galaxies Honig et al. (2010) constrained the num- 
ber of clouds along equatorial rays and their distribution 
using high-angular resolution ground-based mid-IR spec- 



troscopic observations, while fixing the other torus model 
parameters. They found however, no differences in the 
number and distribution of clouds between type 1 and 
type 2 AGN. 

In this paper we expand on our previous work by com- 
bining photometric SEDs with high angular resolution 
(~ 0.3 — 0.4'^) mid-IR spectroscopic observations of a 
sample of 13 nearby Seyfert galaxies. For the first time 
in this work we fit these data for a sizeable sample of 
Seyferts to put tighter constraints on torus model pa- 
rameters. This paper is organized as follows. Section 2 
presents the new observations and the data compiled 
from the literature. Section 3 describes the clumpy dusty 
torus models and the modelling technique. A discussion 
of the fits and inferred torus parameters is presented in 
Section 4. Section 5 discusses the properties of the torus 
and AGN. Finally our conclusions are given in Section 6. 
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Fig. 1. Sketch of the CLUMPY models of Nenkova et al. 
(2008a, b). The radial torus thickness Y is defined as the ratio 
between outer radius (Ro) and the dust sublimation radius (Rd)- 
All the clouds have the same ry , and crtorus characterizes the width 
of the angular distribution of clouds. The number of cloud encoun- 
ters is a function of the viewing angle i, the width of the angular 
distribution of the clouds crtorus, and the mean number of clouds 
along radial equatorial rays A^o (see Equation 2). 



2. OBSERVATIONS 

2.1. The sample 

Our sample is composed of 13 nearby Seyfert galax- 
ies at a mean distance of 38Mpc (median of 31Mpc). 
We selected the galaxies to have high (< 0.8^') angu- 
lar resolution near-IR (1 — 5 /im) imaging, and mid-IR 
(8 — 20 /im) imaging and spectroscopy in the literature 
(Sections 2.3 and 2.4). Although the sample is not com- 
plete and may not be representative, the high- angular 
resolution IR data ensure that we can isolate the IR 
emission arising from the torus surrounding the AGN. 
Currently such high angular resolution imaging and spec- 
troscopy can only be attained from the ground using 8- 
lOm-class telescopes (typically ~ 0.3 — 0.4^' at 10 /im). 
The properties of this sample relevant to this work, and 
the corresponding references are summarized in Table 1. 
The sample includes type 1 and type 2 Seyfert galaxies, 
as well as Seyfert 2 galaxies with a broad component de- 
tected in the near-IR. The bolometric luminosities of the 
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TABLE 1 
Properties of the sample. 



Galaxy 


z 


d 


Type 


Ref 


Nh (X-ray) 


Ref 


€)cone 


Ref 


Sl2/xm 


Ref 


log Lboi 


Ref 






Mpc 






lO^^cm ^ 




deg 




pc 




erg s ^ 




Circinus 


0.001448 


4 


Sy2 


Al 


430. 


Bl 


80-90 


CI 


0.4,2.0 


Dl 


43.6 


El 


IC 4329A 


0.016054 


65 


Syl.2 


A2 


0.61 


B2 


— 


— 


<10.8 


D2 


45.0 


E2 


IC 5063 


0.011348 


46 


Sy2 


A3 


21.78 


B2 


60 


C2 


~ 




44.3-44.7 


E2 


MCG —5-23-16 


0.008486 


34 


by 2 (broad Pap) 


A4 


1.6 


B3 






2.8 


D2 


44.4 


E3 


NGC 1068 


0.003793 


15 


Sy2 


A5 


> 1000. 


B4 


40, 65 


C3 


0.5 X 1.4, 3x4 


D3 


45.0 


E4 


NGC 2110 


0.007789 


31 


Sy2 (broad Br7) 


A6 


2.84 


B2 


30 


C4 






43.8-43.9 


E2 


NGC 3227 


0.003859 


17 


Syl.5 


A5 


1.74 


B2 


60-70 


C5 






43.2-43.5 


E2 


NGC 4151 


0.003319 


13 


Syl.5 


A5 


6.9 


B5 


67, 75 


C3 


2 


D4 


43.7 


E4 


NGC 5506 


0.006181 


25 


NLSyl (broad Pa^) 


A4, A7 


2.78 


B2 


90-100 


C6 




D2 


44.1-44.3 


E2 


NGC 7172 


0.008683 


35 


Sy2 


A2 


8.19 


B2 










43.7-43.8 


E2 


NGC 7469 


0.016317 


66 


Syl 


A5 


0.05 


B6 






10.5 


D2 


45.0-45.1 


E2 


NGC 7582 


0.005254 


21 


Sy2 (broad Br7) 


A6 


5.,^ 100. 


B7 


86, 120 


C3 




D2 


43.3 


E2 


NGC 7674 


0.028924 


118 


Sy2 (broad Pa/3) 


A8 


> 1000. 


B8 




C2 






- 45 


E5 



Notes. — Distances are for ifo = 75kms -"^Mpc ^ and for the nearby objects are taken from RA09 and RAIL References. — Al. Oliva 
et al. (1994). A2: Veron-Cetty & Veron (2006). A3: Colina, Sparks, & Macchetto (1991). A4: Blanco et al. (1990). A5: Osterbrock 
& Martel (1993) and references therein. A6: Reunanen et al. (2003). A7: Nagar et al. (2002). A8: Ruiz et al. (1994). Bl: Matt et al. 
(1999). B2: Vasudevan et al. (2010). B3: Perola et al. (2002). B4: Risaliti, Maiohno, & Salvati (1999). B5: Beckmann et al. (2005). 
B6: Guainazzi et al. (1994). B7: Bianchi et al. (2009). B8: Matt et al. (2000). CI: Maiohno et al. (2000). C2: Schmitt et al. (2003). 
C3: Wilson & Tsvetanov (1994). C4: Pogge (1989). C5: MundeU et al. (1995). C6: Wilson et al. (1985), Maiohno et al. (1994). Dl: 
Tristram et al. (2007). D2: Tristram et al. (2009). D3: Raban et al. (2009). D4: Burtscher et al. (2009). El: Moorwood et al. (1996). 
E2: Vasudevan et al. (2010). E3: From the 2 — lOkeV flux of Weaver & Reynolds (1998) and applying a bolometric correction of 20. E4: 
Woo & Urry (2002) and references therein. E5: Estimated from the scattered 2 — lOkeV luminosity by Malaguti et al. (1998). 



AGN, LboL are taken from the works of Woo & Urry 
(2002) and Vasudevan et al. (2010), or estimated from 
the hard X-ray luminosities using a typical bolometric 
correction of 20 (Elvis et al. 1994). The AGN bolomet- 
ric luminosities of our sample span almost two orders of 
magnitude, from ^ 10^^ to 10"^^ ergs~^. In terms of their 
X-ray neutral hydrogen column density, A^h (X-ray), the 
sample contains both Compton thin objects and Comp- 
ton thick galaxies. We finally list in Table 1 two observa- 
tional properties related to the torus, namely the open- 
ing angle of the ionization cones Ocone, and the torus size 
'5i2/im ( FWHM) as derived from the modelling of mid-IR 
inter ferometric observations. These two parameters will 
be compared with the fitted torus model parameters in 
Section 5.1. 

As explained in the Introduction, in this paper in ad- 
dition to the IR photometric SEDs we will fit ground- 
based mid-IR (~ 8 — 13 \im) spectroscopy, including the 
9.7 /im silicate feature. Shi et al. (2006) found that the 
9.7 /im silicate features measured from S'pitzer jW^^ spec- 
troscopy of a large sample of AGN vary from emission 
to absorption with increasing neutral hydrogen column 
densities. These authors interpreted this result with a 
scenario where the obscuring material is located in two 
different physical scales (0.1 — 10 pc disk and a large 
scale disk extending up to 100 pc) with the dust dis- 
tributed in clouds of different properties. Although the 
mid-IR ground-based spectroscopy in this work (see Sec- 
tions 2.2 and 2.4) probes much smaller physical scales 
than those probed by the 5pi^2;er/IRS data, it is of in- 
terest to place our sample of galaxies in the context 
of larger samples. In Figure 2 we show the apparent 
strength of the 9.7 /im silicate feature (S9.7) as a func- 
tion of the X-ray hydrogen column density for Seyfert 
Is, Seyfert 2s, PG quasars, and 2MASS quasars adapted 
from the work of Shi et al. (2006). The strengths of 
the silicate feature in this figure are mostly measured 



from 5pi^2;er/IRS spectroscop}{3 and thus correspond to 
large physical sizes (at the mean distance of our sam- 
ple - 2kpc). For NGC 4151, NGC 7469 and NGC 3227 
5*9 7 were also measured from S'pitzerjW^^ spectroscopy 
(Thompson et al. 2009), and for NGC 7582, NGC 7674 
and IC 5063 from ground-based data (Honig et al. 2010, 
and Section 2.4). From this figure it is clear that our 
relatively small sample probes well the observed ranges 
for Seyfert Is and Seyfert 2s. We will come back to this 
issue in Section 5. 

2.2. New observations 

A^-band spectroscopy of NGC 4151 was obtained with 
Michelle (Glasse et al. 1997) on the Gemini North tele- 
scope on 20070318 (Program ID GN-2006B-Q-18). The 
standard mid-IR chop-nod technique was employed for 
the observations, and the chop and nod distances were 
15'^ The LowN grating and 2-pixel (0.36^^) slit were used, 
giving R ~ 200. The slit was orientated at 60 degrees 
E of N, along the extended mid-IR emission discovered 
by Radomski et al. (2003). Two blocks of spectroscopy 
were obtained, each of 450 sec on-source (in the guided 
chop beam) , and the galaxy nucleus was reacquired in be- 
tween. Telluric standard stars were observed before and 
after the NGC 4151 spectroscopy, with the same obser- 
vational setup. Michelle data files contain planes consist- 
ing of the difference image for each chopped pair for each 
nod. The chopped pairs were examined for anomalously 
high background or electronic noise, but none needed to 
be excluded from the final data set. The chopped dif- 
ference images were then combined until a single file was 
obtained for each spectroscopy block, and this file for the 
two blocks then averaged together. The resulting spec- 
trum was extracted in a 2-pixel aperture, wavelength- 
calibrated using telluric lines, divided by the standard 
star, and multiplied by a black-body spectrum. 

The typical extraction apertures of the Shi et al. (2006) IRS 
short-low (SL) spectra were 6pixels or 10.8^^ 
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Fig. 2. — Apparent strength of the 9.7//m sihcate feature vs. the 
X-ray hydrogen column density, adapted from the work of Shi et 
al. (2006). For the strength of the silicate feature, positive values 
mean the feature is observed in emission and negative values in 
absorption. Star-like symbols are PG quasars, asterisks are 2MASS 
quasars, circles are Seyfert Is, and squares are Seyfert 2s. The 
filled symbols indicate the galaxies in our sample. The triangles 
are those Seyfert 2s in our sample with broad lines detected in the 
near-IR (see Table 1). The X-ray column densities are from Shi et 
al. (2006) except for the galaxies in our sample, which are taken 
from the references given in Table 1. The strength of the silicate 
features are measured from Spitzer/IRS data and ground-based 
data (see Section 2.1). 



Q-band (~ 20 jam) spectroscopy of NGC 1068 was 
obtained with Michelle on 20090923 (Program ID GN- 
2009B-Q-58), using the same chop- nod technique. The 
3-pixel (0.54^^) slit and lowQ grating were used, for 
R ~ 100, with the slit orientated along the ionisation 
cones (20 degrees E of N). The galaxy nucleus was ob- 
served for a total of 1800 sec on-source, and telluric stan- 
dard stars were observed before and after NGC 1068. 
The data reduction procedure was similar to that used 
for the NGC 4151 A/'-band spectroscopy, including ex- 
traction of the spectrum in a 2-pixel (0.36^^) aperture. 
However, cancellation of the strong and variable Q-band 
telluric lines was poor in the initial, ratioed spectrum. 
This was improved by adding a sloping baseline to the 
standard star spectra before dividingF^. The effect on 
the spectral slope of the resulting NGC 1068 spectrum 
is small. The flux calibration was done with the telluric 
standard stars. The resulting flux densities are probably 
affected by slit losses, and therefore in Section 2.4 we 
recalibrate the spectra with photometric observations. 
Finally the Q-band spectrum used for the modelling was 
rebinned using a 3 pixel box. 

2.3. Published unresolved nuclear fluxes 

i.e., fitting a continuum, multiplying by a few per cent, adding 
that back to the star spectrum, and dividing by that instead. 



As discussed at length in our previous papers (Alonso- 
Herrero et al. 2001, 2003; RA09 and RAll) high angular 
resolution observations are required to isolate the emis- 
sion associated with the torus, and with the direct view 
of the AGN in type Is as well. At the distances of our 
galaxies and the current angular resolutions of the near 
and mid-IR imaging and spectroscopic observations the 
torus emission appears unresolved. In the near-IR up 
to A 2 /im, extended stellar emission arising in the 
host galaxy contaminates and even dominates the nu- 
clear fluxes of type 2 Seyferts (Alonso-Herrero, Ward, & 
Kotilainen 1996; Videla et al. 2011) and is not negligible 
in type 1 Seyferts (Kotilainen et al. 1992). At longer 
wavelengths (A > 3 /im) contamination by stellar pho- 
tospheric emission is greatly reduced. However, any ex- 
tended nuclear emission not directly related to the dusty 
torus, such as dusty clouds in the NLR and the coronal 
line region (e.g.. Bock et al. 2000; Alloin et al. 2000; 
Radomski et al. 2003; Packham et al. 2005; Mason et 
al. 2006; Roche et al. 2006; Reunanen et al. 2010) 
and/or dust heated by young massive stars (Siebenmor- 
gen, Kriigel, & Spoon 2004; Alonso-Herrero et al. 2006; 
Mason et al. 2007; Reunanen et al. 2010), needs to be 
removed. 

To isolate as much as possible the emission of the torus 
(and the AGN when seen directly) we compiled high- 
angular resolution near-IR and mid-IR fluxes from the 
literature with estimates of the unresolved emission when 
available. These unresolved fluxes are the result of re- 
moving, using various methods, the underlying near-IR 
stellar emission and the mid-IR emission produced by 
star formation. The compiled photometry includes near- 
IR ground-based and iJS'T/NICMOS observations (with 
the NIC2 camera, angular resolutions 0.15 — 0.2^') and 
ground-based mid-IR measurements with angular reso- 
lutions of < 0.8'^ The only exception are the L-band 
measurements of Ward et al. (1987) for type 1 Seyferts 
(Table 1), but we use them as upper limits. All the mid- 
IR photometric points have angular resolutions in the 
range 0.3 — 0.5^' to match the resolution and slit widths 
of the mid-IR spectroscopic data (see Section 2.4). When 
possible, for a given galaxy we tried to match the angular 
resolution of the photometric points to make sure we are 
modelling similar physical scales. Table 2 lists for each 
galaxy in our sample the wavelengths and references of 
the photometric data used to construct their SEDs. Fi- 
nally, we used the S'pitzerjW^^ 30 /im continuum fluxes of 
Deo et al. (2009) for those galaxies in our sample without 
nuclear star formation as upper limits in our fits. 

Based on discussions in the papers listed in Table 2 
and comparisons between different works we use the fol- 
lowing errors for our analysis. For near-IR ground-based 
data, except for the NACO data (see below) of Prieto 
et al. (2010), we use for the J-band 30%, for the R 
and K bands 25%, and for the L-band 20%. These in- 
clude the photometric error, the background subtraction 
uncertainty and the uncertainty from estimating the un- 
resolved flux. For the last one we note that the stellar 
emission contribution within a given aperture decreases 
with increasing wavelength and thus the stellar contri- 
bution has a minimum in the L-band (see e.g., Koti- 
lainen et al. 1992 and Alonso-Herrero et al. 1996). 
The M-band fluxes are always considered as upper limits 
because estimating the unresolved component was not 
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TABLE 2 

Wavelength coverage of the observations and references 



Galaxy 


near-IR SED 


Ref. 


mid-IR SED 


Ref 


Spectroscopy 


Ref 


Circinus 


JHKK'L'M' 


Al 


8.7, 18.3 /im 


Bl 


AT-band MIDI 0.60'' 


CI 


IC 4329A 


HK 


A2 


8.99, 11.88 /im 


B2 


AT-band VISIR 0.75'' 


C2 




L 


A3 










IC 5063 


H 


A4 


8.7, 18.3 /im 


Bl 


A^-band T-ReCS 0.67" 


C3 




K 


A5 


10.5, 11.3, 11.9, 12.3 /im 


B2 






MCG -5-23-16 


JHKL'M 


A2 


8.59, 8.99, 11.88 /im 


B2 


A^-band VISIR 0.75" 


C2 








18.72 /im 


B3 






NGC 1068 


JHK 


A2 


8.8, 18.3 /im 


B4 


AT-band Michelle 0.36" 


C4 




LM 


A6 






Q-band Michelle 0.54" 


This work 


NGC 2110 


L' 


A7 


11.2 /im 


B5 


A-band Michelle 0.36" 


C5 








8.6, 9.0, 11.9/im 


B2 












20 /im 


B6 






NGC 3227 


HK 


A8 


11.3 /im 


Bl 


AT-band VISIR 0.75" 


C2 




L 


A3 


8.99, 11.88 /im 


B2 






NGC 4151 


JHK 


A8 


10.8, 18.2 /im 


Bl 


A-band Michelle 0.36" 


This work 
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References. Al: Prieto et al. (2010). A2: Alonso-Herrero et al. (2001). A3: Ward et al. (1987), the photometry is used as an upper limit. 
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et al. (2007). 



possible. For the NACO AO observations of Prieto et 
al. (2010), which were measured through 0.1-0.2^' aper- 
tures, the smaller contamination by stellar emission when 
compared to natural seeing observations results in lower 
uncertainties. We therefore use 20% in J, and 15% in 
the HKLM NACO bands. Finally, for the NICMOS 
data, which have very stable photometric calibration and 
PSFs, we use 20% in the J-band, and 10 - 20% in the 
H and K bands (unless otherwise specified in the cor- 
responding references). The estimated NICMOS uncer- 
tainties are based on the comparison of the unresolved 
fluxes reported for the same galaxies in different works 
(Alonso-Herrero et al. 2001; Quillen et al. 2001; Gal- 
limore & Matthews 2003, and Kishimoto et al. 2007). 
For the N and Q band measurements we use 15% and 
25% errors, respectively to account for the photometric 
calibration and unresolved component uncertainties (see 
details in RA09). 

We finally address the issue of possible variability of the 
near-IR fluxes and the simultaneity of the SEDs. Nine 
out of the 13 galaxies in our sample have been reported to 
show variability in the near-IR: NGC 7674 (Quillen et al. 
2000), NGC 1068, IC 4329A, NGC 2110, MCG -5-23- 
16, NGC 5506, and NGC 7469 (Glass 2004), NGC 3227 
(Suganuma et al. 2006), and NGC 4151 (Koshida et al. 
2009). There are no reports on mid-IR variability of our 
sources. The typical variability in the near-IR is 40% 
on average, with amplitude variations around the me- 
dian of between 10 — 30% for our sample. For all these 
galaxies the 1 — 2.2 /im data, when available, were taken 
simultaneously. The LM fluxes in some cases were not 



simultaneous with the shorter wavelength fluxes. How- 
ever, in most cases the LM fluxes are taken as upper 
limits as there was no estimate of the unresolved emis- 
sion. Glass (2004) showed that the L-band variability is 
typically less than 20% around the median flux, which 
is within the photometric and unresolved emission un- 
certainties. We can therefore assume that variability in 
the near-IR does not affect the compiled SED within the 
above discussed uncertainties of the unresolved measure- 
ments. 

2.4. Published ground-based mid-IR spectroscopy 

The published mid-IR spectroscopy used in this work 
(see Table 2) was obtained with four different instru- 
ments: 

• (i) The Thermal-Region Camera Spectrograph (T- 
ReCS; Telesco et al. 1998) on the Gemini-South 
Telescope. The T-ReCS data were obtained with 
the low-resolution grating, which provides a spec- 
tral resolution of ^ 100. The slit widths were 
0.67'' (Young et al. 2007) and 0.36'' (Roche et al. 
2006, 2007). 

• (ii) The Michelle instrument (Glasse et al. 1997) on 
the Gemini-North telescope. The Michelle obser- 
vations (Mason et al. 2006; 2009 and Section 2.2) 
were obtained with a 0.36" slit width and a spec- 
tral resolution oi R ^ 200 in the A^-band, and a 
0.54" slit and R - 100 in the Q-band. 

• (iii) VISIR, the mid-IR imager and spectrograph 
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Fig. 3. — The black lines are the original flux-calibrated Michelle 
and Q-band spectra of NGC 1068 (see Mason et al. 2006 and 
Section 2.1, respectively). The large and small square symbols 
are the photometry of Tomono et al. (2001) for a circular 0.4^^- 
diameter aperture on the original data, and for a 0.29^^ x 0.15^^ 
aperture on the deconvolved data, respectively. The green lines are 
the and Q-band spectra scaled to the 0.4^^ photometric points at 
8.7 and 18.5 /xm photometric points. After the scaling, the Q-band 
spectrum was also rebinned (see text for details). 



0.4'' photometric points. Finally, for all the A/'-band 
spectra in our sample we added in quadrature the intrin- 
sic error of the spectra and the 15% error of the imag- 
ing data (see Section 2.3) used for scaling the data, for 
each spectroscopic point. For the Q-band spectrum of 
NGC 1068 the associated error for scaling the spectrum 
is 25%. These errors were added for the fitting. 

3. MODELING OF THE DATA 

3.1. Clumpy Torus Models 

In this work we use an interpolated version (see Sec- 
tion 3.3) of the clumpy dusty torus models of Nenkova et 
al. (2002, 2008a, b) including the corrections for the pre- 
viously erroneous AGN scaling factor (see the erratum 
by Nenkova et al. 2010). The CLUMPY models are de- 
scribed by the six parameters listed in Table 3 that we 
explain in the following (see also Figure 1). An AGN with 
a bolometric luminosity Lboi(AGN) is surrounded by a 
torus of dusty clouds and all the clouds have the same 
optical depth ry, which is defined in the optical F-band. 
The torus clouds are located between the inner radius of 
the torus and the outer radius of the torus Rq^ and 
the torus radial thickness is defined as Y = Rq/R^. The 
inner radius of the torus is set by the dust sublimation 
temperature (Tgub ^ 1500 K), 



mounted on the 8.2m UT3 telescope at the 
ESO/Paranal observatory in Chile. The VISIR ob- 
servations (Honig et al. 2010) were obtained in 
low-spectral resolution mode {R ^ 300) with a slit 
width of 0.75''. 

• (iv) MIDI (Leinert et al. 2003), the mid-IR in- 
terferometer at the VLTI at the ESO/Paranal ob- 
servatory. The MIDI observations (Tristram et al. 
2007) were taken with a 0.6" width slit and R ~ 30. 
The data used here are a "total fiux spectrum", 
where 4 spectra are obtained for this, two for each 
telescope of the interferometer and 2 for each "win- 
dow" in MIDI. The spectrum was extracted in a 6 
pixel wide mask or 0.516", and calibrated individ- 
ually for the windows and telescopes. 

Because the VISIR ground-based mid-IR spectra have 
slightly higher spectral resolution than the other spectra, 
we rebinned them to contain approximately 150 spectral 
points. The MIDI, T-ReCS and Michelle spectra contain 
between 100 and 200 spectral points. We did not at- 
tempt to remove any emission lines (e.g., [Siv]10.51 fim) 
or broad features (e.g., PAH features) in the spectra. 

The final step was to scale the spectra to the photo- 
metric points. This was mostly necessary for the spectra 
observed with the narrowest slits (see Table 2) because 
of slit losses. As an example, in Figure 3 we show the 
original and Q-band spectra (Mason et al. 2006 and 
Section 2.1, respectively) of NGC 1068 together with the 
mid-IR photometric points of Tomono et al. (2001) for a 
circular aperture of 0.4"-diameter, and the photometry 
through a rectangular 0.29" x 0.15" aperture after decon- 
volution of the data. To match the angular resolutions 
of the imaging and spectroscopic data, the NGC 1068 TV 
and Q-band spectra were scaled to their corresponding 



The angular distribution of the clouds is assumed to 
have a smooth boundary and it is described as a Gaus- 
sian with a width parameter atorus- The radial distribu- 
tion is a declining power law with index q (ex r~^). The 
mean number of clouds along radial equatorial rays is 
Nq. The number of clouds along the line of sight (LOS) 
at a viewing angle i (measured from the polar direction, 
see Figure 1) is, 

NLOsii) = No e-(^°-^)'/^'orus. (2) 

In a clumpy dust distribution the classification of an 
object as a type 1 or type 2 AGN is not truly a mat- 
ter of the viewing angle but of the probability for direct 
view of the AGN (see Figure 1, and also Elitzur 2008). 
The probability for an AGN-produced photon to escape 
through the torus along a viewing angle i when all the 
clouds are optically thick (ry > 1) is. 

Peso ^ 6^-'^^°^^ (3) 

In the models the radiative transfer equations are 
solved for each clump and thus the solutions depend 
mainly on the location of each clump within the torus, 
its optical depth, and the chosen dust composition. We 
adopt a dust extinction profile corresponding to the 
OHMc dust (i.e., the standard cold oxygen-rich ISM dust 
of Ossenkopf et al. 1992). The total torus emission is 
calculated by integrating the source function of the total 
number of clumps convolved with the radiation propaga- 
tion probability along the torus (Nenkova et al. 2002). 
For Type 1 and intermediate type Seyferts where there 
is an unobscured view of the AGN, it is also possible to 
include its contribution to the resulting SED. The AGN 
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TABLE 3 

Axial ratios, foreground extinction measurements and 
strength of the silicate feature 
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Notes. — The axial ratios are the ratio of to minor to major axis of 
the host galaxies and are taken from de Vaucouleurs et al. (1991), 
except for that of NGC 7172 that is from Jarrett et al. (2003). 
S9.7 are the observed strengths of the silicate feature, with posi- 
tive numbers indicating that the feature is in emission, and neg- 
ative numbers the feature is in absorption. For the Roche et al. 
(2007) galaxies, the two measurements correspond to fits to the 
feature done using two different silicate grain profiles. References. 
Al: Maiohno et al. (2000). A2: Storchi-Bergmann et al. (1999). 
A3: Heisler & Robertis (1999). A4: Veilleux et al. (1997). A5: 
Goodrich et al. (1994). A6: Winge et al. (2000). A7: Riffel, 
Rodriguez-Ardila, & Pastoriza (2006). Bl: Roche et al. (2007). 
B2: Honig et al. (2010). B3: Mason et al. (2006). B4: Thompson 
et al. (2009). 

continuum emission in these models is characterized with 
a piecewise power law distribution (see Nenkova et al. 
2008a for details). 

We finally note that there is some evidence of the pres- 
ence of an extra hot dust component, originating very 
close to the AGN. Reverberation mapping (see discus- 
sion by Kishimoto et al. 2007) and near-IR interfer- 
ometric observations (e.g., Kishimoto et al. 2009 and 
references therein) of Seyfert galaxies have shown that 
the subhmation radii appear to be smaller than expected 
from Equation 1 (similar to the Barvainis 1987 relation). 
However in this work we take the simplest approach for 
fitting the data assuming that all the near-IR emission is 
originating in the classical torus, with no additional hot 
dust components. 

3.2. Foreground dust components 

In this section we examine the possibility of having con- 
tributions to the observed IR SEDs and mid-IR spectra 
of Seyfert galaxies from dust emission and absorption un- 
related to the AGN. There are several pieces of evidence 
we can look into, namely the inclination of the galax- 
ies, and the presence of nuclear dust lanes and spatially 
resolved variations of the silicate feature. 

Approximately half of the galaxies in our sample have 
minor-to-major axis ratios of b/a < 0.5 (inclined host 
galaxies, see Table 4), and as suggested by Deo et al. 
(2009), it is likely that dust in the host galaxy disk can 
contribute significantly to the observed silicate absorp- 
tion and the long-wavelength continuum. The only ex- 
ception appears to be IC 4329A, which is a highly in- 
cHned system, but does not show silicate absorption. Deo 
et al. (2009) concluded for this galaxy that it is likely 
that our fine of sight does not intersect any dense clouds 



TABLE 4 

Parameters of the CLUMPY Torus Models 



Parameter Symbol Interval 



Torus radial thickness Y [5, 30] 

Torus angular width crtorus [15°, 70°] 

Number of clouds along an equatorial ray Nq [1, 15] 

Index of the radial density profile q [0, 3] 

Viewing angle i [0°, 90°] 

Optical depth per single cloud ry [5, 150] 



Notes. — Torus radial thickness: Y = Ro/Rd where Ro is the outer 
radius and R^ is the inner radius (Equation 1). The cloud distri- 
bution between R^ and Ro is parametrized as r~^. 

in the host galaxy. 

Using color maps IVIartini et al. (2003) showed that in 
general the dust in the nuclear regions of Seyfert galax- 
ies is on physical scales not associated with those of the 
dusty torus. It is then likely that this dust is located in 
front of the nucleus (see also Regan & IMulchaey 1999) 
and arises from the galactic ISIM. In our sample, Circinus, 
IC 5063, NGC 5506, NGC 7582, NGC 2110, NGC 7172, 
and possibly NGC 7674 show dust features in the nuclear 
(central ~ 1 — 2^') regions (see e.g., Quillen et al. 1999; 
Maiolino et al. 2000; Martini et al. 2003). 

Some of the galaxies in our sample show variations of 
the 9.7 /im silicate feature on scales of 1 — 2^' indicat- 
ing the presence of extended dust components (Mason 
et al. 2006; Roche et al. 2006, 2007: Colling, Roche, 
& Mason 2009). Similarly, the spatially-resolved mid- 
IR polarimetric observations of NGC 1068 (Packham et 
al. 2007) can be explained with a geometrically, and op- 
tically thick torus surrounded by a larger, more diffuse 
structure, associated with the dusty central regions of 
the host galaxy. 

Roche et al. (2007) showed that NGC 5506 and 
NGC 7172, both with prominent nuclear dust lanes (see 
e.g., Malkan et al. 1998), have similar absorbing columns 
as derived from X-ray observations and the 9.7 /im fea- 
tures, for the latter just by using a foreground dust screen 
model. This suggests that some of the extinction mea- 
sured from the 9.7 /im silicate feature might arise in the 
galactic ISM (see also the discussion in Section 2.1, and 
Figure 2). In other words, in these two galaxies if the sil- 
icate feature were to come only from the torus, it would 
be shallower (at it is filled in by emission from warm dust 
also in the torus). These findings would also be consis- 
tent with the work of Guainazzi, Matt, & Perola (2005) 
who showed that the presence of dust lanes on scales of 
100 pc in Compton-thin Seyfert 2s is correlated with the 
X-ray obscuration, mostly in the X-ray column density 
range ~ 10^^ — 10^^ cm~^. This was interpreted by these 
authors as due to the larger covering fraction of the gas 
in the dust lanes, rather than the pc-scale dusty torus. 

From a theoretical point of view, clumpy dusty models 
cannot produce very deep silicate features {Sg,7 < —1) 
(see Nenkova et al. 2008b; Elitzur 2008; Sirocky et al. 
2008), while observations show that many Seyfert galax- 
ies have relatively deep silicate features (see Shi et al. 
2006; Hao et al. 2007; Deo et al. 2009). In our sample, 
this includes Circinus, NGC 5506, NGC 7582, NGC 7172 
(see Table 4 and Figure 2). By analogy with those 
ULIRGs optically classified as type-2 AGN and with deep 




A. Alonso-Herrero et al. 

5 

NGC3227 



100 




100 



Rest-frame wavelength (/xm) 



Rest-frame wavelength (fim) 




100 



Rest-frame wavelength (/xm) 



Rest-frame wavelength (fim) 



Fig. 4. — Model fits for pure Seyfert 1 galaxies with torus + AGN emission. The filled dots are the SED photometric data and the black 
line is the mid-IR spectrum. The fits of the torus emission to the IR SED data points alone are shown in orange, while the fits to the IR 
SED and spectroscopy are shown in blue. The solid blue and orange lines are the models described by the combination of parameters that 
maximizes their probability (MAP) distributions, while the dashed lines correspond to the model computed with the median value of the 
probability distribution of each parameter. The shaded areas indicate the range of models compatible with a 68% confidence interval. The 
inset shows in detail the fit to the A^-band spectroscopy. We did not restrict any of the torus model parameters. No foreground extinction 
was used for the fits. 



silicate features modelled by Sirocky et al. (2008), the 
deep silicate features of some Seyfert 2 galaxies in our 
sample may be explained by additional obscuration by 
cold foreground dust. 

For simplicity we will consider the extended dust com- 
ponent as a purely foreground medium producing only 
absorption. We note however, that high angular resolu- 
tion mid-IR spectroscopic observations of the ionization 
cone dust of NGC 1068 show evidence of components of 
both absorption and emission. We use the IR extinc- 
tion curve of Chiar & Tielens (2006) of the local ISM 
in the wavelength range ~ 1 — 35 /im, which includes 
the two silicate features at 9.7 and 18 /im. For extinc- 
tions {Ay (frg) > 5 mag) the effects of foreground extinc- 
tion cannot be ignored, especially in the spectral region 
around the 9.7 /xm silicate feature. 

Finally, it is important to note the possible degeneracy 
between AGN type and the effects of foreground extinc- 



tion. Using smooth density torus models Alonso-Herrero 
et al. (2003) demonstrated that the reddened near-IR 
SED from a nearly face-on (polar view) torus (underly- 
ing Seyfert 1) looks similar to the intrinsic near-IR SED 
at a viewing angle closer to a true Seyfert 2 (equato- 
rial view) given sufficient levels of foreground extinction 
(Ay (frg) > 5 mag). Given this degeneracy in our fitting 
process we choose not to leave the foreground extinc- 
tion as a free parameter (see Section 3.3), but instead we 
use pubhshed values of the extinction as estimated from 
near-IR measurements. Table 4 gives the values of the 
foreground extinction used in this work, only for galax- 
ies with evidence of extended dust components and with 
Ay (frg) > 5 mag. The foreground extinction is applied 
to the torus emission, as well as to the AGN emission 
that is included for galaxies with a direct view of the 
BLR. 
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Rest-frame wavelength (yLtm) Rest-frame wavelength (fim) 

Fig. 5. — Model fits for pure Seyfert 2 galaxies with torus emission alone. Symbols and lines are as in Figure 4. Upper left panel: 
Circinus. We fixed the foreground extinction to Ay (fig) = 9 mag (Maiolino et al. 2000). We restricted the viewing angle to the range 
i = 60 — 90°. Upper right panel: IC 5063. We did not restrict any of the torus model parameters. The foreground extinction was fixed to 
A^(frg) = 7 mag based on results by Heisler de Robertis (1999). Lower left panel: NGC 1068. We restricted the viewing angle to the 
range i = 60 — 90° (see Section 3.3). For NGC 1068 the second inset shows in detail the fit to the Q-band spectroscopy. Lower right panel: 
NGC 7172. Only the SED+spectroscopy fit is shown due to the limited number of photometric points in the near-IR. We did not restrict 
any of the torus model parameters. The foreground extinction was fixed to Ay (frg) = 40 mag. 



3.3. Modeling Technique 

The CLUMPY database now contains 1.2 x 10^ mod- 
els, which are calculated for a fine grid of model param- 
eters. To fit the data we take a Bayesian approach to 
dealing with the inherent degeneracy of the torus model 
parameters (see RA09 and RAll). In this work we use 
an updated version of the BayesClumpy code developed 
by Asensio Ramos & Ramos Almeida (2009). This tool 
uses a Bayesian inference to allow as much information as 
possible to be extracted from the observations. We refer 
the reader to Asensio Ramos & Ramos Almeida (2009) 
for details on the interpolation methods and algorithms 
used by BayesClumpy. However, it is important to note 
that in using a Bayesian approach we do not make use 
of the original set of models of Nenkova et al. (2008a, b), 
but rather of an interpolated version of them (see figures 
3 and 4 in Asensio Ramos & Ramos Almeida 2009). The 



fineness of the grid of clumpy models makes interpolation 
an appropriate methodology for our studies. 

The new version of BayesClumpy allows, in addition 
to fitting photometric points, for the possibility of fit- 
ting spectra. In the case of photometric observations, 
BayesClumpy simulates the transmission curve of the 
corresponding filter on the model spectra. For the spec- 
troscopic observations, the full information, flux and as- 
sociated uncertainty, is used. When both photomet- 
ric SEDs and spectroscopy are fitted together, we use 
the data uncertainties discussed in Section 2 as weights 
for the fitting routine. Under the assumption that the 
torus models are valid simultaneously for photometric 
and spectroscopic data and that the noise in all observed 
points is mutually uncorrelated, the Bayesian approach 
is insensitive to the fact that there are many more spec- 
troscopic points than photometric ones. The reason is 
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Fig. 6. — (a) Model fits for Seyfert 2 galaxies with broad lines detected in the near-IR. Upper left panel: MCG —5-23-16. Model fits are 
for torus emission alone. All symbols and lines are as in Figure 4. The viewing angle was restricted to a Gaussian distribution centered 
at i = 55° with a 10° width (see Section 3.3). The foreground extinction was fixed to Ay(frg) = 7 mag based on results by Veilleux 
et al. (1997). Upper right panel: NGC 2110. Model fits are for torus + AGN emission. The viewing angle was restricted to a Gaussian 
distribution centered at i = 40° with a 10° width (see Section 3.3). Only the SED+spectroscopy fit is shown due to the limited number 
of photometric points in the near-IR. The foreground extinction was fixed to Ay(frg) = 5 mag (Storchi-Bergman et al. 1999). Lower left 
panel: NGC 5506. Model fits are for only for torus emission (see Section 4.2). The viewing angle was restricted to a Gaussian distribution 
centered at i = 40° with a 10° width (see Section 3.3). The foreground extinction was fixed to A\/(frg) = 11 mag (Goodrich et al. 
1994). Lower right panel: NGC 7582. Model fits are for torus emission alone. We did not restrict any of the torus model parameters. The 
foreground extinction was fixed to Ay(frg) = 13 mag (Winge et al. 2000). 



that, apart from the possible regularization that can be 
included in the prior distributions, we sample the full 
posterior distribution and obtain marginalized posterior 
for each parameter. The marginalization procedure takes 
into account all possible values of the parameters produc- 
ing good fits. This avoids the possible overfitting of the 
SED which is produced by neglecting the few photomet- 
ric points. 

The prior distributions for the model parameters are 
assumed to be truncated uniform distributions in the 
ranges given in Table 3. We note that in the most up- 
to-date version of the CLUMPY models after correct- 
ing for the erroneous AGN scaling factor (Nenkova et 
al. 2010), the optical depth of the individual clouds only 
goes up to 150, instead of ry = 200 of the older mod- 



els. The only prior information we use in this work is 
the viewing angle. For those galaxies in our sample with 
H2O maser detections: Circinus (Greenhill et al. 2003) 
and NGC 1068 (Greenhih et al. 1996) we restricted the 
viewing angles i to values in the range 60 — 90°, that 
is, close to equatorial views through the torus. We can 
also use the accretion disk viewing angles deduced from 
X-ray observations of the FeK a line as an additional 
constraint, if we assume that the accretion disk and the 
torus are coplanar. We found estimates for three galaxies 
in our sample, for MCG —5-23-16 is z ~ 53° (Reeves et 
al. 2007), and for NGC 5506 and NGC 2110 is i - 40° 
(Guainazzi et al. 2010 and Weaver & Reynolds 1998). 
For these three galaxies we assumed Gaussian distribu- 
tions for this parameter with a width of 10 °. 
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TABLE 5 

Fitted torus model parameters from SED 
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Fig. 6. (b) Continued. NGC 7674. Model fits are for torus 
emission alone. We did not restrict any of the torus model param- 
eters. The foreground extinction was fixed to Av(frg) = 5 mag. 



In addition to the the six torus model parameters, 
there are two extra parameters that can be fitted or fixed. 
The first parameter accounts for the vertical displace- 
ment needed to match the fluxes of a given model to an 
observed SED/spectrum. This vertical shift, which we 
allow to vary freely, scales with the AGN bolometric lu- 
minosity (see Nenkova et al. 2008b) and will be discussed 
in Sections 4.1 and 4.2. The second parameter is the fore- 
ground extinction (see Section 3.2 and Table 4) due to 
the host galaxy, which is different from that produced by 
the torus along the LOS (see RA09 and RAll). 

For the modelhng of the SED and spectroscopy of type 
1 Seyferts the AGN contribution needs to be added to the 
torus emission (see Nenkova et al. 2008a for the assumed 
shape). In principle, the same should be done for those 
Seyfert 2s in our sample with broad lines detected in the 



near-IR, as this means we have a direct view of the BLR. 
However most of these Seyfert 2s in our sample are also 
heavily affected by extinction (see Table 4). Kishimoto 
et al. (2007) demonstrated for type 1 Seyferts that most 
of the unresolved emission at 2.2 fim is produced by hot 
dust emission from the inner walls of the torus, and that 
there is a very small contribution from the big blue bump 
emission (i.e., AGN emission) at this wavelength. Then 
there is the question for AGN with broad lines and fore- 
ground extinctions Ay ^ 5 mag of whether we are seeing 
AGN emission in the range 1 — 2 /im or not. We thus de- 
cided for these galaxies to fit the data both including and 
not including the AGN emission. This will be discussed 
in more detail in Section 4.2. 

4. RESULTS FROM FITS TO 
SED+SPECTROSCOPY 

The result of the fitting process are the marginal poste- 
rior distributions for the six free parameters that describe 
the CLUMPY models plus the vertical shift. These are 
the probability distributions of each parameter, which 
are represented as histograms. As explained in Sec- 
tion 3.1, except for the viewing angle of five galaxies, 
for the rest of the parameters we use uniform priors. If 
the observational data contain sufficient information for 
the fit, then the resulting probability distributions of the 
fitted torus model parameters will clearly differ from uni- 
form distributions. In those cases the probabilities either 
show trends or are centered at certain values within the 
considered intervals. 

For each galaxy we translate the fitted torus model 
parameters into two model spectra. The first one corre- 
sponds to the maximum- a-posteriori (MAP) values that 
represent the best fit to the data. The second is the 
model produced with the median value of the proba- 
bility distribution of each parameter, which is charac- 
teristic of the observed SED+spectroscopy data. Fig- 
ures 4 to 6 show for Seyfert Is, Seyfert 2s, and Seyfert 
2s with broad lines in the near-IR, respectively, these fits 
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Fig. 7. — Marginal posterior distributions of the free parameters that describe the CLUMPY models resulting from fitting the data for 
IC 4329A The SED alone fits are shown in orange and the SED+spectroscopy fits are shown in blue. The solid and dashed lines indicate 
the MAP and median values of the distributions, and the shaded areas are the ±lcr values. We did not use a foreground extinction. The 
marginal posterior distributions for the rest of the sample are shown in Appendix A. 



to the SED+spectroscopy. Figure 7 shows the marginal 
posterior distributions of the six torus model parame- 
ters for IC 4329A. The marginal posterior distributions 
for the rest of the sample are shown in Appendix A. 
In Appendix B we also show for Circinus, NGC 1068, 
NGC 4151, and IC 4329A the two-dimensional poste- 
rior distributions for all combinations of the torus model 
parameters. This kind of two-dimensional distributions 
can be used to check for possible correlations and degen- 
eracies between different torus model parameters for a 
given galaxy. Table 5 summarizes the statistics for the 
fitted torus model parameters. In Sections 4.1 to 4.3 we 
discuss new constraints on the torus model parameters 
when fitting the SED+spectroscopy data together, and 
in Section 4.4 we assess the improvements obtained by 
adding mid-IR spectroscopy to the SED data. 

4.1. Fits to pure Seyfert Is and Seyfert 2s 

In this section we discuss the fits to the pure Seyfert 
Is and 2^, while the fits to those Seyfert 2s with broad 
lines detected in the near-IR (including NGC 5506) are 
discussed separately in Section 4.2. As can be seen from 
Figures 4 and 5 (in blue), the CLUMPY torus models 
provide very good simultaneous fits to the photometric 
SED and spectroscopy data of Seyfert Is and 2s, in par- 
ticular for those galaxies with low host galaxy foreground 
extinction. 

It is worth noting that all the Seyfert Is except 

In this work pure type 1 AGN have broad lines detected in 
the optical (i.e., include types 1.5, 1.8 and 1.9), while pure type 
2s do not have broad lines detected either in the optical or in the 
near-IR. 



IC 4329A show a slight excess of emission in the near-IR, 
above the median torus+AGN model fits, which might 
be attributed to hot dust. Mor et al. (2009) included, 
apart from the torus emission, two extra components, 
hot dust and NLR emission, to fit the Spitzer/lUS spec- 
tra of PC quasars. These extra components provided the 
additional flux need in the near-IR for their sample. It is 
not clear however if such components are needed in our 
fits, because the unresolved emission we used for our fits 
probes typically a few tens of par sees, while the Mor et 
al. (2009) data cover physical sizes on scales of a few kpc. 
Finally, the near-IR photometric points of our Seyfert Is 
are well within the ila confidence regions of the fitted 
models. 

The only galaxy for which we could not get a good 
fit to both the SED photometric points and the mid-IR 
spectroscopy was the Seyfert 2 galaxy NGC 7172. This 
galaxy is highly inclined (see Table 4), has prominent 
dust lanes, and thus probably suffers from high extinc- 
tion in the host galaxy (see e.g. Roche et al. 2007). 
Since we could not find estimates for the host galaxy ex- 
tinction, we used the foreground extinction derived from 
the X-ray column density and the standard Galactic ra- 
tio N^/Av = 1.910^^cm"2mag"^ (Bohlin, Savage, & 
Drake 1978). As can be seen from Figure 5, while the 
model provides a reasonable fit to the 9.7 jiia silicate fea- 
ture, the near-IR photometric points are well above the 
model. One possibility is that there is contamination by 
extended dust components in this galaxy, and the un- 
resolved flux estimates from the K and L-band ground- 
based data are upper limits. However, the most likely 
explanation is that the clumpy torus emission + cold 
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foreground dust screen are not the appropriate model 
for this galaxy. Rather, a spherically symmetric smooth 
model may be more appropriate for deeply embedded 
objects, such as NGC 7172 (see Levenson et al. 2007). 

As can be seen from the marginal posterior distribu- 
tions in Figure 7 and Figures Al to A7 (blue lines and 
blue shaded regions), the majority of the fitted torus 
model parameters of the pure Seyfert Is and Seyfert 2s 
are well constrained. In particular, the width of the an- 
gular distribution atorus, the radial torus thickness F, 
and the viewing angle z, which were also relatively well 
constrained from fits to the SED alone (RA09, RAll, 
and Section 4.4), present narrow probability distribu- 
tions. As found by RAll there is no clear relation be- 
tween the derived viewing angle i and the classification 
of the galaxy into a type 1 or a type 2. In other words, 
not all type Is are viewed at relatively low inclination 
angles and not all type 2 are seen at directions closer to 
the equatorial plane of the torus. However, as we shall 
see in Section 5.2, the relevant quantity for a galaxy to 
be classified as a type 1 or a type 2 is the probability for 
an AGN produced photon to escape unabsorbed or not. 

As can be seen from Table 5, our results confirm that 
CLUMPY models of tori with relatively small radial 
thicknesses, produce good fits with no need for very large 
tori. The fitted values are within the assumed range of 

Y = 5 — 30, and the data for most galaxies can be fit- 
ted with tori with radial thickness Y ^ 10 — 15. In 
Section 5.1 we will compare in more detail our derived 
physical sizes of the torus (in parsecs) with those derived 
from the modeling of mid-IR interferometry. 

The index of the power law distribution q controls the 
placement of clouds between the inner, hotter parts of 
the torus and the outer cooler regions (Nenkova et al. 
2008b), and thus plays a role in the shape of the IR SED 
(see also Honig & Kishimoto 2010). The effects of chang- 
ing q on the CLUMPY models are easier to see in type 
1 objects (low values of i). For steep radial distributions 
(clouds concentrated near the inner radius of the torus) 
the near and mid-IR SEDs become redder (see figure 9 of 
Nenkova et al. 2008b), although including the AGN con- 
tribution always makes the model SEDs flatter. We find 
that the fitted values of q for our Seyfert Is are relatively 
small. However the reverse for Seyfert 2s (seen at high 
inclinations) is not necessarily true. We also note that 
the sensitivity of the SED to q for very small values of 

Y is highly reduced. This is because for such small tori 
the SED does not change much whether the clouds are 
distributed along the whole extent of the torus {q = 0) 
or highly concentrated in its inner part (large q). This 
lack of sensitivity could result in posterior distributions 
that depend on the quality of the interpolation technique 
used by BayesClumpy. Consequently, whenever the in- 
ferred Y is very small (e.g., NGC 1068), the inferred 
value of q should be treated with caution. 

For a larger sample and using fits of the SEDs alone 
RAll found statistically significant differences of q be- 
tween type Is and type 2s. RAll found that the SEDs 
of type Is were fitted with small values of while those 
of type 2s required larger q. Honig et al. (2010), on 
the other hand, found the q parameter (|a|, in their no- 
tation) to be in the range — 1.5 for both type Is and 
2s. The other torus model parameters, Nq and ry, de- 
rived from the SED+spectroscopy data will be discussed 
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Fig. 8. — Comparison between the AGN bolometric luminosities 
(median values) as derived from the SED+spectroscopy fitting and 
those from the literature. The solid line is the 1:1 relation, not a fit. 
Filled symbols as in Figure 2. For those galaxies with broad lines 
and high extinctions we show the inferred AGN bolometric lumi- 
nosities from fits with and without the AGN component. The for- 
mer are indicated as "(+AGN)". Note that the derived bolometric 
luminosities for NGC 2110 with and without the AGN component 
are the same. 



in the context of the fits of the spectral region around 
the 9.7 /im silicate feature (Section 4.3). For a detailed 
discussion of the fits to individual sources and compar- 
ison with clumpy torus models fits in the literature, we 
refer the reader to Appendix B. 

As explained in Section 3.3, the shift applied to scale 
the CLUMPY models to the data is directly related to 
the bolometric luminosity of the AGN. In this work we 
chose to leave the shift as a free parameter, but other 
works (e.g., Mor et al. 2009) used it as a constraint for 
the fits. Figure 8 shows a comparison between the fitted 
AGN bolometric luminosities (^[^^^^^^(AGN)) and AGN 
bolometric luminosities from the literature. The latter 
were computed with different methods, including apply- 
ing bolometric corrections and modelling of the SEDs 
(see Table 1 and Section 2.1). The typical uncertain- 
ties of the AGN bolometric luminosities from scaling the 
CLUMPY models are 0.1 — 0.2dex. It is clear from this 
figure that the agreement between the bolometric lumi- 
nosities is good for the majority of pure Seyfert Is and 
Seyfert 2s, with most of our estimates within ~ 0.4 dex 
of the 1:1 relation. The only exception is the Seyfert 
1 NGC 7469. The fitted AGN bolometric luminosity for 
this galaxy is below two independent literature estimates 
(e.g.. Woo & Urry; Vasudevan et al. 2010). There is 
strong nuclear (~ 0.2^' ~ 65 pc) star formation in this 
galaxy, but it only contributes a small fraction of the K- 
band luminosity within 0.2^' (see e.g., Davies et al. 2007 
and references therein). However, the ground-based mid- 
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IR nuclear spectrum of this galaxy (Honig et al. 2010) 
show faint 11.3 /im PAH emission, and the PAH emis- 
sion becomes very prominent in the circumnuclear re- 
gions (Roche et al. 1991). It is then likely that the mid- 
IR nuclear fluxes and spectroscopy of NGC 7469 contain 
a contribution from star formation. A similar situation 
may be the case for NGC 3227. 

4.2. Fits to Seyfert 2s with near-IR broad lines 

Four galaxies in our sample are classified as Seyfert 2s, 
but there are also reports in the literature of detections 
of broad lines in the near-IR (see Table 1 for references). 
Additionally, NGC 5506 has different spectral classifi- 
cations in the literature and there is some controversy 
about whether it has broad components of the near-IR 
emission lines (Blanco et al. 1990; Ruiz et al. 1994: 
Goodrich et al. 1994; Veilleux et al. 1997). Nagar et al. 
(2002) clearly detected BLR emission in the near-IR and 
classified NGC 5506 as a narrow line Seyfert 1 (NLSyl). 
Given these facts we discuss NGC 5506 in this section. 

As the near-IR data suggest, for this kind of galaxies we 
may have a direct view of the BLR, and thus in principle 
we should include the AGN component when modelling 
the data. However, some of these galaxies also tend to 
suffer from relatively high values of foreground extinction 
(Table 4). For these Seyferts we performed the fits with 
and without the AGN components. In this section we use 
other observational properties of these galaxies to deter- 
mine whether the AGN component should be included 
or not for fitting their SED+spectroscopy data. 

In Section 4.1 we showed that the model fits for pure 
Seyfert Is and 2s provided good constraints to the AGN 
bolometric luminosity. We can then use the compari- 
son between the fitted AGN luminosity and those taken 
from the literature to determine if we need to add the 
AGN component. Figure 8 suggests that for NGC 5506, 
NGC 7582, NGC 7674, and MCG -5-23-16 we do not 
need to add the AGN component. For NGC 2110 we 
obtained comparable bolometric luminosities (and other 
torus model parameters) with and without the AGN 
component. This is probably due to the lack of pho- 
tometric points in the near-IR for this galaxy. Finally 
we note that when we included the AGN component to 
fit the SED -h spectroscopy data of NGC 5506 we could 
not fit the data with the viewing angles inferred from 
X-ray data (Section 3.2). This is in agreement with the 
luminosity comparison. 

The fits of the Seyfert 2s with broad lines detected 
in the near-IR are shown in Figure 6, while the marginal 
posterior distributions are displayed in Appendix A (Fig- 
ures A8 to A12). For the three galaxies with moder- 
ate silicate features (MCG -5-23-16, NGC 2110, and 
NGC 7674, see Table 4) we are able to reproduce reason- 
ably well the photometric SED points and the mid-IR 
spectroscopy, and the torus model parameters are mostly 
well constrained. As for the pure Seyfert 2 NGC 7172, 
the two galaxies with relatively deep silicate features in 
absorption (NGC 5506 and NGC 7582), it is possible to 
do a simultaneous fit of the photometry and mid-IR spec- 
troscopy, although it is not a good fit. Imanishi (2000) 
measured an unusual ratio between the depths of the 
3.4 /im carbonaceous dust absorption and the 9.7 /im sil- 
icate dust in NGC 5506. He suggested that the obscura- 



tion towards the nucleus of NGC 5506 could be ascribed 
mostly to dust in this host galaxy (>100-pc scale) and 
not to the obscuring torus. This may also reflect differ- 
ences in grain populations (see Roche et al. 2007). The 
X-ray emission of NGC 7582 indicates the presence of at 
least two absorbers of columns densities of ~ 10^^ cm~^ 
and 4 — 5 X 10^^ cm~^, with the latter probably associated 
with large scale obscuration (Bianchi et al. 2009). 

4.3. Fits to the silicate feature 

The insets of the upper panels of Figures 4 to 6 show 
in detail the fits in the spectral region around the 9.7 /im 
silicate feature. For NGC 1068 we also show in Figure 5 
the fit to the Q-band spectrum. It is clear that the inter- 
polated version of the Nenkova et al. (2008a, b) clumpy 
torus models produce excellent fits to the A^-band silicate 
feature, and the ~ 20 /im spectrum of NGC 1068. In par- 
ticular for galaxies with low host galaxy extinctions and 
moderate silicate strengths {89,7 > —1), the fits to the 
silicate feature spectral region are also compatible with 
those to the SED photometric points. This suggests that 
the data mostly correspond to the torus emission. These 
galaxies are ah the pure Seyfert Is (IC 4329A, NGC 3227, 
NGC 4151 and NGC 7469), and NGC 1068, IC 5063, 
MCG -5-23-16, and NGC 7674. 

We now discuss the fits of those galaxies with deep sil- 
icate features. For the Circinus galaxy, which has a rela- 
tively deep silicate feature {89,7 = —1.8/ — 2.4, Roche et 
al. 2007 and Table 4), CLUMPY torus models with the 
parameters given in Table 5 plus a foreground absorb- 
ing screen with Ay {fig) = 9 mag reproduce reasonably 
well the silicate feature and the photometric SED, except 
for the 18.3 /im flux. For the other galaxies (NGC 5506, 
NGC 7172 and NGC 7582), although the fits to the fea- 
ture are qualitatively good, the models do not reproduce 
at the same time all the photometric points. One possi- 
bility for these three galaxies is that there is an impor- 
tant contamination from extended dust structures, while 
in Circinus, the closest galaxy in our sample (d = 4 Mpc) 
this contamination is minimized. 

Nenkova et al. (2008b) emphasized that the clumpy 
torus distributions produce more elaborate patterns of 
the 9.7 /im and the 18 /im silicate features (see also Honig 
& Kishimoto 2010), while for smooth density models the 
silicate feature is always in emission for face-on views and 
in absorption for edge-on views. Besides, the CLUMPY 
torus models never produce very deep silicate features, in 
contrast with smooth density models (e.g.. Pier & Krolik 
1993; Efstathiou & Rowan- Robinson 1995). When fitted 
together with the photometric SEDs we expect the sili- 
cate feature to be mostly sensitive to the optical depth 
of the clouds ry , and the average number of clouds along 
radial equatorial rays Nq. As an illustration of the com- 
plicated behavior of the silicate feature in terms of the 
torus model parameters, we can see that galaxies of the 
same type and with similar apparent depths of the 9.7 /im 
silicate feature (e.g., type Is IC 4329 A and NGC 7469, 
and type 2s MCG-5-23-16 and IC 5063, see Honig et 
al. 2010) have different fitted values of Nq and ry (see 
Table 5, and Appendix A). 

A general result for our sample of Seyfert galaxies is 
that the average number of clouds along radial equatorial 
rays is never very low, and is in the range A^o = 8 — 15. 
This is in good agreement with the results of Nenkova et 
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TABLE 6 

AGN AND TORUS MODEL PROPERTIES DERIVED FROM THE FITS. 
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Notes. — L™^^^^^(AGN) and rtorus were derived from the median values of the marginal posterior distributions. 



al. (2008b). However, Honig et al. (2010) inferred fewer 
clouds along equatorial rays from their fits to the mid- 
IR spectroscopy of Seyfert galaxies, but we note these 
authors fixed the value of atorus? among other param- 
eters (see Appendix B for a more detailed discussion). 
We also note here, that the values of the optical depth 
of the clouds ry in the corrected version of the models 
(see Nenkova et al. 2010) only go up to 150, whereas 
the older version were up to 200. In some cases (e.g., 
NGC 3227, NGC 7469) it appears as if to compensate 
for limited values of ry, the fits are achieved with more 
clouds along radial equatorial rays, close to the maximum 
value allowed by the models of A^o = 15. 

While for most Seyferts in our sample CLUMPY torus 
models with a typical value of A'o = 12 produced good 
fits to the data, the Spitzer/lH^ spectra of PC quasars 
were well fitted with torus models containing a mean 
value of A^o = 5 (Mor et al. 2009). This is weh under- 
stood because most PG quasars show the 9.7 /im silicate 
feature in emission (Shi et al. 2006; Hao et al. 2007), and 
the Nenkova et al. clumpy torus models with atorus ^ 30° 
and A'o > 10 almost always produce the feature in ab- 
sorption, for all viewing angles (see also discussions by 
Nenkova et al. 2008b and Nikutta et al. 2009). 

The optical depths of our sample of Seyferts show a 
broad distribution, ry ~ 50 — 150, with no obvious de- 
pendence on other torus model parameters or AGN type 
(see also Mor et al. 2009). This result may be under- 
stood because the CLUMPY model near and mid-IR 
SEDs and the strength of silicate silicate feature vary 
only slightly with varying i for a given set of A'o and 
cTtorus, at TV > 100 (see Nenkova et al. 2008b). The IR 
data of the PG quasars studied by Mor et al. (2009) were 
fitted with TV < 100 (Mor et al. 2009). 

4.4. Comparison with fits of SEDs alone 

In Figures 4-6 we also compare the fits done with 
SED+spectroscopy (blue lines and blue shaded regions) 
and SEDs alone (orange lines and orange shaded re- 
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Fig. 9. — Width of the torus crtorus vs. the complementary to the 
viewing angle /3 = 90 — i. We plot the median and the ±lcr values of 
the probability distributions of the fitted values. The dashed lines 
are lines of constant AGN escape probabilities Peso (Equation 3) 
for an average number of clouds along an radial equatorial rays of 
A^o = 12, which is representative of our sample (see Section 4.3). 
Filled symbols as in Figure 2. 



gions), except in cases of sparse SEDs. Some of our 
galaxies are in common with the works of Ramos Almeida 
et al. (2009, 2011). We have redone the SED alone fits 
for those galaxies to use consistent errors with our work. 
In general we find that the main difference is that fits 
to SEDs alone tend to infer broader probability distribu- 
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tions of these torus model parameters (see Appendix A). 
In particular, for galaxies with low foreground extinc- 
tions the fits with SEDs alone when compared to the 
fits with SED+spectroscopy produce compatible values 
(within the ±lcr values of the probability distributions) 
of CTtorus, and i. 

Most disagreements in the fitted parameters tend to 
occur for A^o, and ry. Some of these differences 
are explicable. For instance, Nq and q may trade off 
to yield an approximately constant number of clouds 
in the inner region. Also, the optical depth of the in- 
dividual clouds Ty and the average number of clouds 
along an equatorial ray A^o have a strong influence on 
the depth of the 9.7 /im silicate feature and the shape 
of the ~ 8 — 13 /im emission. Note also that the SED 
alone fits of some pure type 1 Seyferts always tend to 
produce the 9.7 /im silicate feature in emission (see Fig- 
ure 4), while the observations show an almost featureless 
spectrum (IC 4329A, NGC 3227) or the feature slightly 
in emission (NGC 4151). For the Seyfert 2s with low host 
galaxy foreground extinction there is a better agreement 
for these two parameters between fits with and without 
mid-IR spectroscopy. 

The other important fact to note is that including 
mid-IR spectroscopy helps assess whether there is an 
important contribution from extended dust structures 
causing absorption (and possibly emission too). This 
is clearly the case for Circinus, NGC 5506, NGC 7172 
and NGC 7582, which have relatively deep silicate fea- 
tures. Indeed, the fits to the SEDs alone of Seyfert 2s 
of RA09, which were done without foreground extinc- 
tion, only predicted moderately deep silicate features for 
these galaxies. However, we reiterate that it is not al- 
ways possible to produced good simultaneous fits to the 
SEDs and spectroscopy of Seyferts with deep silicate fea- 
tures even when the host galaxy foreground extinction is 
included. It is possible that for the Seyfert galaxies with 
the deepest silicate features a clumpy medium in a torus- 
like configuration may not be appropriate to explain the 
observations (see e.g., Levenson et al. 2007). 

5. PROPERTIES OF THE TORUS AND AGN 
5.1. Torus size and Angular width 

In this section we discuss the two torus model param- 
eters that can be compared with observations, the torus 
size and the angular width. 

For each AGN, using the Y parameter and the bolo- 
metric luminosity of the system derived from the fits, and 
Equation 1 we computed the physical radius rtorus of the 
torus. We used the median value of the fitted torus radial 
thickness Y and AGN bolometric luminosity. As can be 
seen from Table 6 the torus radii derived from our fits to 
the SED+spectroscopy are between r torus ^ 

1 and 6pc. 

Our fitted physical sizes are consistent with the mid-IR 
interferometric results that the torus is relatively com- 
pact (Tristram et al. 2007, 2009; Burstscher et al. 2009; 
Raban et al. 2009). 

Before we compare the torus radii with the 12 /im sizes 
inferred from mid-IR interferometric observations, it is 
necessary to discuss some caveats. First, the modelling 
of the 12 /im interferometric observations requires struc- 
tures on different physical scales, as demonstrated for 
Circinus and NGC 1068 (see Tristram et al. 2007, 2009; 



Raban et al. 2009, and Table 1). However, the FWHM 
sizes of the cooler, more extended component could still 
be compared with model predictions. Second, Nenkova 
et al. (2008b) showed that for tori with F = 10 - 30 and 
various viewing angles, at 12 /im ~ 60 — 75% of the total 
flux is enclosed within angular radii of 6 ^ 2 — 56 d- Here 
6>d is the angular size at the inner radius of the torus. 
This is because the 12 /im emission traces warm dust, 
with T > 200 K, and is rather insensitive to cooler ma- 
terial further from the nuclei. Finally the compact dusty 
torus is expected to blend into the galaxy disk (e.g., Shi 
et al. 2006; Packham et al. 2007; Nenkova et al. 2008b), 
and thus deriving the true size of the torus might not be 
straightforward. 

With these considerations in mind, we can compare the 
sizes derived from the modelling of the 12 /im interfero- 
metric data may not be appropriate. However, values of 
Y larger than the 12 /im sizes from interferometric obser- 
vations are a necessary, although not sufficient, condition 
for consistency between our modelling and the interfero- 
metric results. There are estimates of the 12 /im physical 
size of the torus through mid-IR interferometry for six 
galaxies in our sample (see Table 1 for sizes and corre- 
sponding references). We find a relatively good agree- 
ment, within a factor of two, between the fitted torus 
sizeq3 and the sizes inferred from mid-IR interferometry. 
The galaxy with the largest discrepancy is NGC 7469. 
Tristram et al. (2009) inferred 5i2/im = 10.5 pc, whereas 
our fitted torus size is about 5pc. We note however, 
that the interferometric observations of NGC 7469 were 
affected by large uncertainties in the determination of 
the total flux that may have compromised the derived 
physical size of the torus of this galaxy. 

The torus model width of the angular distribution of 
the clouds, (Jtorus, is related to the opening angle of the 
ionization cones, Ocone- The latter can be measured from 
emission line images, usually [O ill] A5007 and/or Ha. Ta- 
ble 1 gives the measured opening angles (measured from 
the polar direction) for our sample and corresponding 
references. It is important to note that this compari- 
son is not straightforward because the relation between 
the true opening angle of the cone and the measured 
opening angle from line emission depends on the gas 
distribution. For thin gas disks, Mulchaey, Wilson, & 
Tsvetanov (1996) simulations showed that cone-shaped 
emission can be produced for most viewing angles and 
that the observed angles are less than the true opening 
angle. We can try to make this comparison for galax- 
ies with the observed widest ionization cones, for which 
we can set more meaningful upper limits to the widths 
of the torus. For the two galaxies with Bcone ^ 100° 
(NGC 5506 and NGC 7582, Table 1) the fitted atoms are 
compatible with the approximate limits set by the obser- 
vations ((Jobs < 40 — 50°). Note however, the difficulties 
encountered for fitting the observations of these galaxies 
(Section 4.2). 

5.2. AGN Escape Probability 

One of most salient properties of the clumpy torus 
models is that there is always a finite probability for an 

For comparison with the 12 /im sizes (FWHM) from mid-IR 
interferometric observations, the rtorus values given in Table 6 need 
to be multiplied by two. 
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Fig. 10. — AGN bolometric luminosities vs. the torus geometrical covering factor /2 (left panel) and AGN bolometric luminosities vs. 
the angular width of the cloud distribution of the torus crtorus (right panel). The quantities for our sample of Seyferts were derived from 
the SED+spectroscopy fits (filled symbols, as in Figure 2). The plotted uncertainties for /2 and crtorus are the ±lcr values of the posterior 
distributions of the fitted parameters. Also plotted are PG quasars (star-like symbols), for which Mor et al. (2009) derived the torus /2 and 
ctorus using the Nenkova et al. (2008b) models. The bolometric luminosities of the PG quasars are derived from the optical luminosities 
(see Mor et al. 2009 for details). 



obscured view of the AGN, irrespective of the viewing 
angle (Nenkova et al. 2008a, b). While in the simplest 
version of the unified model type Is would be observed 
in directions closer to the pole of the torus (low values of 
i) than type 2s, in a clumpy medium this is not necessar- 
ily so. Indeed, it is possible to have a type 1 viewed at a 
relatively high i provided that the width of the angular 
size of the torus is small (for a given A/'o), as is the case 
for NGC 4151 and IC 4329A. Note that two-dimensional 
posterior distributions of these two galaxies show that 
the fitted values of i and (Jtorus are inversely correlated 
(Figures B3 and B4 in the Appendix.) Seyfert 2s can also 
be viewed from any angle, as long as there is at least one 
cloud along the line of sight. The probability of encoun- 
tering a cloud is higher for increasing values of crtorus, 
A^o, and i. The fact that not all type 2s in our sample re- 
quire very broad angular distributions (large crtorus, e.g., 
NGC 1068 in Figure 9), suggests that Nq and i must be 
high in those cases. In Figure 9 the error bars reflect the 
±lcr uncertainty around the median of the probability 
distributions of the fitted parameters. 

A more relevant quantity for an AGN to be classified 
as a type 1 or a type 2 is the probability for an AGN- 
produced photon to escape unabsorbed (Equation 3). 
Since the escape probability is a highly non-linear func- 
tion of A/'o, i, and crtorus, we took full advantage of our 
Bayesian approach. We thus generated the posterior dis- 
tributions for Pesc from the posterior distributions of the 
relevant fitted parameters given in Table 5. As can be 
seen from the values in Table 6, all pure type Is and type 
2s have respectively, relatively high (~ 12 — 44%) and low 
(< 0.1%) AGN escape probabilities, as expected. The 
type 2 Seyferts with broad-lines detected in the near-IR 



show a range of escape probabilities. We note that these 
probabilities do not include the effects of foreground ex- 
tinction. 

5.3. Covering Factors and a Receding Torus? 

Observational (e.g., Simpson 2005) and theoretical 
(e.g., Honig & Beckert 2007) arguments provide evidence 
of the so-called receding torus (Lawrence 1991). In this 
scenario the higher fraction of type 1 AGN at high AGN 
luminosities is explained in terms of a decreasing cover- 
ing factor of the torus. Furthermore, Honig & Beckert 
(2007) and Nenkova et al. (2008b) argued that the de- 
creasing fraction of type 2s at high luminosities not only 
arises from decreasing crtorus, but also from decreasing 
A/'o, or both. 

In this work we calculated two different covering fac- 
tors using our Bayesian approach. The first one is the 
geometrical covering factor /2 , which does not depend on 
the viewing angle, and is derived by integrating the AGN 
escape probability (Equation 3) over all angles as put for- 
ward by Mor et al. (2009). If we define /3 = 7r/2-z, then: 

.7r/2 

/2 = 1- / Pesc{P)cOs{p)dp. (4) 

A related quantity is the the ratio between the torus 
luminosity integrated over the entire wavelength range 
covered by the torus models (^ 0.2 — 700 /im) and the 
inferred AGN bolometric luminosity: 

i:model(AGN)- 

This ratio is the apparent covering factor (see Mor et 
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al. 2009), because it depends on the viewing angle i (as 
well as on atoms and A^o), and can be understood as a 
reprocessing efficiency. As we did for Pesc we generated 
the posterior distributions for /2 and f{i) from the poste- 
rior distributions of the relevant fitted parameters given 
in Table 5. The median values and the ±lcr values of the 
distributions of the geometrical covering factor and re- 
processing efficiencies of our sample of Seyferts are given 
in Table 6. 

The galaxies in our sample show a range of reprocessing 
efficiencies, although in general there is no clear depen- 
dency with any torus model parameters. In contrast, the 
torus geometrical covering factor /2 appears to depend on 
the AGN bolometric luminosity, as can be seen from Fig- 
ure 10 (left panel). We also included in this comparison 
high luminosity AGN. These are PG quasars, for which 
Mor et al. (2009) computed the torus /2 from Nq and 
c^torus as inferred from fits to Spitzer/IRS spectroscopy. 
We note that these authors also used the Nenkova et 
al. (2008a, b) torus models, but added extra components 
to fit their data (see Section 4.1). Mor et al. (2009) 
calculated the PG quasars AGN bolometric luminosities 
using the observed optical luminosities (see Mor et al. 
2009 for more details). For our sample of Seyfert galax- 
ies the uncertainties in /2 take into account the ila 
values around the median of the probability distribution 
of 

<^torus- When combining the PG quasars of Mor et 
al. (2009) and our sample of Seyferts, there is a ten- 
dency for the torus covering factor /2 to decrease from 
/2 0.9-1 at low AGN luminosities (lO"^^-"^"^ ergs"^) to 
/2 ~ 0.1 — 0.3 at high luminosities (> 10^^ ergs~^). Such 
tendency for the torus covering factor to decrease with 
the bolometric luminosity of the AGN was already noted 
for PG quasars by Mor et al. (2009) when including an 
additional hot black-body component to fit their data. 

A similar trend may also be present when comparing 
atoms vs. L;^°^^^(AGN), with high luminosity AGN hav- 
ing narrower cloud angular distributions (Figure 10, right 
panel) than low luminosity ones. We also found in Sec- 
tion 4 that on average Seyfert galaxies (with bolomet- 
ric luminosities in the range 10^^~^^ ergs~^) are fitted 
with more clouds along a radial equatorial ray than PG 
quasars (bolometric luminosities > lO^^ergs"^). Given 
these trends, it is not surprising that the torus geomet- 
rical covering factor /2, which depends on both Nq and 
c^torus, tends to be lower in the most luminous objects in 
our sample than in the less luminous ones. 

The decreased torus covering factor at high AGN lumi- 
nosities may provide some support for a receding torus 
at high luminosities. We note however, that in these 
comparisons there are few type 2 objects at high AGN 
luminosities, so we cannot rule out a possible depen- 
dency of atoms with the AGN type. That is, type 2s 
might have wider tori than type Is as suggested by the 
results of RAll. Finally a note of caution. A large frac- 
tion of the low luminosity AGN in our sample with wide 
angular cloud distributions are in highly inclined galax- 
ies. Therefore, we cannot exclude the possibility that 
contamination from extended dust features not directly 
related to the dusty torus has resulted in large fitted val- 
ues of atoms- Finally, as we explained in Section 2.1, our 
sample is not complete, and may not be representative, 
and thus this result needs to be explored for a complete 



sample. 

6. SUMMARY AND CONCLUSIONS 

This is the third paper in a series performing detailed 
fits to the unresolved IR emission of AGN using the 
clumpy torus models of Nenkova et al. (2008b). In 
the first two papers of the series (RA09 and RAll) we 
fitted only the photometric SEDs, while in this paper 
we also included high angular resolution (0.3 — 0.4'') 
ground-based mid-IR spectroscopy. The sample in this 
work is composed of 13 nearby (at a median distance 
of 31 Mpc) Seyfert galaxies with bolometric luminosities 
Lboi(AGN) - 10"^^ - lO'^^ergs-^ The sample contains 
pure type 1 and type 2 Seyferts, as well as Seyfert 2s 
with broad-lines detected in the near-IR. The galaxies 
also span a range in the observed strength of the 9.7 fim 
silicate feature, going from mild emission to deep absorp- 
tion, and in X-ray column densities, from Compton thin 
to Compton thick objects. 

We compiled near and mid-IR imaging and spec- 
troscopy data taken with high-angular resolution (~ 
0.3 — 0.8''). We also presented new ground-based mid-IR 
spectroscopic observations in the TV-band of NGC 4151 
and in the Q-band of NGC 1068. We have used an inter- 
polated version of the clumpy torus models of Nenkova et 
al. (2008a, b and 2010) and a Bayesian approach to fit to- 
gether the photometric SEDs and mid-IR spectroscopy. 
The fits provided the probability distributions for the 
six torus model parameters that describe the CLUMPY 
models (see Table 3), plus the AGN bolometric lumi- 
nosity. We included the effects of foreground extinc- 
tion in the host galaxy for galaxies with evidence of 
extended dust structures and Ay (fig) > 5 mag. When 
compared to fits done using SEDs alone, the fits to the 
SED+spectroscopy data result in smaller uncertainties 
(narrower probability distributions) for the angular dis- 
tribution of the clouds atoms, torus radial thickness F, 
and viewing angle i. 

The interpolated version of the Nenkova et al. 
(2008a, b; 2010) clumpy torus models provided good fits 
to the photometric SEDs and mid-IR spectroscopy of our 
sample, and in particular for those Seyfert galaxies with 
low or moderate foreground extinctions (^^(frg) ^ 5 — 
10 mag). For three AGN with relatively deep 9.7 /im sili- 
cate features and hosted in inclined galaxies with promi- 
nent dust lanes (NGC 5506, NGC 7172, and NGC 7582), 
we were not able to fit simultaneously the photomet- 
ric SEDs and mid-IR spectroscopy. These three galax- 
ies show evidence of extended dust structures on tens 
of parsec scales, not directly related to the dusty torus. 
In contrast, we were able to fit reasonably well the data 
of Circinus. One possibility is that because Circinus is 
the closest galaxy in our sample (distance of 4 Mpc), the 
contamination by extended dusty structures is reduced. 

CLUMPY models with small tori (torus radial thick- 
ness in the range Y ^ 10 — 15) provided adequate fits 
to the data. Combining the modelled AGN bolometric 
luminosities and the values of F, we inferred that the 
physical radii of the tori of Seyfert galaxies are between 
~ 1 pc and 6pc. For the six galaxies in our sample with 
12 /im sizes derived in the literature from the modelling 
of mid-IR inter ferometric observations, we found a rea- 
sonably good agreement. 

The Nenkova et al. (2008b) models were also able 
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to produce good fits to mid-IR high angular resolution 
spectroscopy, in terms of the shape of the silicate fea- 
ture and continuum slope and not only the apparent 
strength of the feature. When combined with the pho- 
tometric SEDs, the mid-IR spectroscopy allowed us to 
constrain the number of clouds along radial equatorial 
rays Nq and the optical depth of the individual clouds 
Ty. We find that the tori of Seyfert galaxies have typi- 
cally A^o ^ 8 — 15, and the optical depth of the clouds are 
in the range ry = 50 — 150. By comparison the tori of 
PG quasars (Lboi(AGN) > lO^^ergs"^) appear to con- 
tain fewer clouds along radial equatorial rays, typically 
TVo ^ 5, with optical depths of ry < 100 (Mor et al. 
2009). 

From the scaling of the models to the data we derived 
the AGN bolometric luminosities, which were in good 
agreement with those derived with other methods in the 
literature. As also found by RAll, the viewing angle i is 
not the only parameter controlling the classification of a 
galaxy into a type 1 or a type 2. For instance, type 2s 
tend to be viewed at directions closer to the equatorial 
plane of the torus, but small values of i are permitted as 
long as there is one cloud along the LOS. This is because 
in clumpy media, there is always a finite probability for 
an AGN photon to escape absorption, and this probabil- 
ity is a function of i, atorus, and A^o- From our fits, we 
found that this probability is relatively high in type Is 
(Pesc ~ 12 — 42%), while it is found to be low in pure 
type 2s (Pesc < 0.1%). 

We finally discussed the possibility of a receding torus. 
We compared the results for our sample of relatively 
low luminosity AGN with those of the PG quasars an- 
alyzed by Mor et al. (2009). Both samples combined 
together span a range of Lboi(AGN) ~ 10^^~^^ ergs~^. 
We found that high luminosity AGN (Lboi(AGN) > 
lO^^ergs"^) tend to have lower torus covering factors 
(/2 ~ 0.1 — 0.3) because they have narrower tori (smaller 
c^torus) and contain fewer clouds along radial equato- 
rial rays A'o- In contrast lower luminosity ones (at 
Lboi(AGN) - lO'^^-'^'^ergs-i) tend to have high cov- 
ering factors (/2 ^ 0.9 — 1). This might explain the 
decreased observed fraction of type 2 AGN at high lu- 
minosities. We note however, that this result may be 
subject to some caveats. In particular most of the AGN 
in our sample at the low end of Lboi(AGN) are hosted 
in highly inclined galaxies, and it is possible that con- 
tamination by extended dust components not related to 



the dusty torus results in larger fitted values of atorus- 
Also, as discussed in Section 2.1, our sample is not com- 
plete, and in particular there are very few type 2s at high 
AGN luminosities. Therefore our results need to be ex- 
plored further, and to this end we have been allocated 
time on the Gran Telescopio Canarias with the mid-IR 
instrument CanariCam (Telesco et al. 2003) for a large 
mid-IR survey of AGN (Levenson et al. 2008). 
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APPENDIX 

MARGINAL POSTERIOR DISTRIBUTIONS 

Here we present the marginal posterior distributions of the fits to the SED+spectroscopy and SED alone for all the 
galaxies in the sample except for IC 4329A, which is shown in Figure 7. 
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Fig. Al. — Marginal posterior distributions of the free parameters that describe the CLUMPY models resulting from fitting the data for 
NGC 3227. The SED alone fits are shown in orange and the SED+spectroscopy fits are shown in blue. The solid and dashed lines indicate 
the MAP and median values of the distributions, and the shaded areas are the zblcr values. We did not use a foreground extinction. 
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Fig. A2.— Same as Figure Al but for NGC 4151. 
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q Ty i [deg] (fixed range) 

Fig. A4. — Same as Figure Al but for Circinus. We fixed the foreground extinction to A^(frg) = 9 mag (Maiolino et al. 2000). We 
restricted the viewing angle to the range i = 60 — 90° . 
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q Ty i [deg] 

Fig. A7. — Same as Figure Al but for NGC 7172. Only the SED+spectroscopy fit is shown due to the hmited number of photometric 
points in the near-IR. We did not restrict any of the torus model parameters. The foreground extinction was fixed to Av(frg) = 40 mag. 
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q Ty i [deg] (restricted) 

Fig. A9. — Same as Figure Al but for NGC 2110. The viewing angle was restricted to a Gaussian distribution centered at i = 40° with 
a 10° width (see Section 3.3). Only the SED+spectroscopy fit is shown due to the limited number of photometric points in the near-IR. 
The foreground extinction was fixed to A)/(frg) = 5 mag (Storchi-Bergman et al. 1999). 



Torus and AGN properties of nearby Seyfert galaxies 



25 




q Ty i [deg] 

Fig. All. — Same as Figure Al but for NGC 7582. We did not restrict any of the torus model parameters. The foreground extinction was 
fixed to A^(frg) = 13 mag (Winge et al. 2000). 
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Fig. A12. — Same as Figure Al but for NGC 7674. Model fits are for torus emission alone. We did not restrict any of the torus model 
parameters. The foreground extinction was fixed to A\/(frg) = 5 mag. 
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APPENDIX 
NOTES ON INDIVIDUAL SOURCES 

In this section we discuss in detail the fits to individual sources and compare them with similar fits done in the liter- 
ature with clumpy dusty torus models. We also show in Figures Bl to B4 the two-dimensional marginal distributions 
of four selected sources: Circinus, NGC 1068, NGC 4151, and IC 4239A. In these figures the contours indicate the 
68% and 95% confidence regions. These two-dimensional distributions can be used to assess whether or not certain 
torus model parameters are correlated. 

Circinus. Circinus is a Seyfert 2 galaxy located at a distance of 4Mpc. The water maser observations of this galaxy 
(Greenhill et al. 2003) trace a warped edge-on disk, thus indicating that the torus is seen at a high inclination. We then 
used the viewing angle as a prior for our fits and restricted it to the range i = 60 — 90 deg. The fit to the photometric 
SED and mid-IR spectroscopy is shown in Figure 5, while the marginal posterior distributions are in Figure A4. As 
can be seen from the latter figure all the torus model parameters are nicely constrained and our estimated value of the 
AGN bolometric luminosity log LJ^^^^^ (AGN) = 43.5 ergs ^ is in excellent agreement with that of Moorwood et al. 
(1996). Figure Bl shows that the fitted torus model parameters of Circinus are not degenerate at the 68% confidence 
level 

Schartmann et al. (2008) fitted the photometric SED and mid-IR spectroscopy of Circinus using their clumpy torus 
models, although due to long computational times they did not explore the entire parameter space of their models. 
They found an appropriate fit for the SED+spectroscopy as well as for the mid-IR interferometric visibilities with 
i = 70 deg, ^open = 65 deg (this is equivalent to our atoms parameter) and a = —0.5 {\a\ is equivalent to our q 
parameter). These are in excellent agreement with our fits using the interpolated version of the Nenkova et al. (2008b) 
models of i = 66^4 deg, atorus = ^^-11 deg, and q = 0.71q 4 (see Table 5). The largest disagreement is in the torus size 
since their modelling favors an outer radius of 30 pc. Our fitted AGN bolometric luminosity sets the dust sublimation 
radius (see Equation 1) of Circinus at 0.07 pc, whereas the fitted radial extent of the torus Y translates into an outer 
radius of rtoms = 2pc (see Table 5). 

NGC 1068. NGC 1068 is a Seyfert 2 galaxy at a distance of 15Mpc. As with Circinus, the water maser 
observations (Greenhill et al. 1996) suggest that the torus is viewed at a high inclination. We used i = 60 — 90 deg as 
a prior for our fits. The fit to the photometric SED and mid-IR spectroscopy is shown in Figure 5. We note that the 
Nenkova et al. (2008a, b) models provide an excellent fit to the photometric SED as well as to the A/'-band and the 
Q-band spectroscopy of this galaxy. The marginal posterior distributions of NGC 1068 are shown in Figure A6, while 
the the two-dimensional marginal distributions are shown in Figure B2. These figures indicate that all the fitted torus 
model parameters for NGC 1068 are nicely constrained. From our fits we derived an AGN bolometric luminosity of 
logLj5^^j^®^(AGN) = 44.8 erg s~^, which is in good agreement with the estimates of Bock et al. (2000) and Vasudevan 
et al. (2010). 

Honig, Prieto, & Beckert (2008) modelled the IR SED and VLTI/MIDI data using the Honig et al. (2006) models. 
Their modeled torus was found to have a line-of-sight inclination of 90 deg with the clouds are distributed according to 
a power law of r~^/^, a half opening angle of the torus of ^ 40 deg. The average number of clouds, which obscure the 
AGN along the line of sight towards the observer, was Nq = 10. Their fitted torus model parameters are in a relatively 
good agreement with our fitted values of i = 88^3 deg, atoms = 26^4 deg, q = 2.21q"3, and Nq = 14^3 (see Table 5). 
Honig et al. (2008) also estimated the diameter of the dust torus to be ~ 1.4 — 2pc, depending on wavelength and 
orientation, whereas our fitted diameter is 4pc. 

Type Is: IC 4329A, NGC 3227, and NGC 7469. These three Seyfert Is are in common with the work of Honig 
et al. (2010). They modelled the mid-IR spectroscopy of a sample of AGN using their newest clumpy torus models 
(Honig & Kishimoto 2010). Honig et al. (2010) showed that the mid-IR spectroscopy provides sensitive constraints 
to the number of clouds along the equatorial direction and the index of power law (a in their notation) of the cloud 
distribution. They fixed other torus parameters such as the opening angle (45 deg), the radial extent to i^out = 150, 
'Tcioud = 50 (see Honig & Kishimoto 2010 for a detailed discussion). We note that their model parameters were also 
chosen to provide good estimates of the VLTI/MIDI observed visibilities (see Honig et al. 2010 for more details). For 
the three Seyfert 1 galaxies in common with our study they find indexes of the power law distribution of |a| ~ — 1 
and Nq = 7.5 — 10, for their assumed typical viewing angles of i = 30 deg. These compare well with our fitted values 
of q ^ 0.2 — 0.9, although we tend to find slightly more clouds along the equatorial line of sight (see Table 5). 

Type 2s: IC 5063, MCG -5-23-16, NGC 2110, and NGC 7674. These four type 2 galaxies in our sample are 
in common with the work of Honig et al. (2010). They modelled the mid-IR spectroscopy of these galaxies with the 
same parameters as the type Is in the previous section but with typical viewing angles of i = 75 deg. Their modelling 
finds indexes of the power law distribution of the clouds |a| ~ — 1.5 and Nq ^ 5. More clouds along the equatorial 
line of sight are allowed in their modelling if the viewing angles for Seyfert 2s are smaller {i ~ 60 deg). Our fitted 
values of q^ except for NGC 2110, are consistent with their highest values of |a|, whereas we find Nq > 10 (median 
values, see Table 5). 
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Fig. Bl. — Circinus. Two-dimensional marginal posterior distributions (joint distributions) for all combinations of the six torus model 
parameters for the SED+spectroscopy fit. The contours indicate the regions at 68% (solid lines) and 95% (dashed lines) confidence levels. 
See Asensio Ramos &; Ramos Almeida (2009) for more details. 
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B3.— As Figure Bl but for NGC 4151. 
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